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ABSTRACT 


This paper addresses the classical problem of printed circuit board routing: that is, the problem of automatic 
routing by a computer other than by brute force that causes the execution time to grow exponentially as a function 
of the complexity. Most of the present solutions are either inexpensive but not efficient and fast, or efficient and fast 
but very costly. Many solutions are proprietary, so not much is written or known about the actual algorithms upon 
which these solutions are based. This paper presents a formal algorithm for routing traces on a printed circuit board. 
The solution presented is very fast and efficient and for the first time speaks to the question eloquently by way of 
symbolic statements. 


NOMENCLATURE 

The following nomenclature do not appear in alphabetical order as each definition is predicated upon 
previous definitions. 

|i the major axis (horizontal or vertical direction) where the sequence {J ; } is associated with |i. 

6 the minor axis (horizontal or vertical direction) where the sequence { M. } is associated with p. 

RM a rectangular matrix of evenly spaced grid points where a node is defined to be an arbitrary grid point 

belonging to RM. 

RM h RM represented as a sequence of horizontal lines composed of grid points. 

RM V RM represented as a sequence of vertical lines composed of grid points (fig. 1). 
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Figure 1. A representation of an internal matrix. 



the number of available nodes from left to right at the node, (p, 1), that belongs to RM h (fig. 1). 

the number of available nodes from bottom to top at the node, (p, 1), that belongs to RM y . 

the i lh designated node on the major axis where sequence { J ; } is monotonic and bounded. That is, 
if (IB, MB) and (IE, ME) are the initial and terminal nodes respectively, then IB = Jj > J 9 > J 3 > 

J 4 J n = IE, or IB = J 1 < J 3 < J 3 < J4 J n = IE. (These inequalities are consistent with the 

horizontal axis as the major axis. Otherwise, the sequence { J i } is bounded by MB and ME, and the 
inequalities will be consistent with the vertical axis as the major axis. Moreover, |Jj — J ; j | repre- 
sents the number of available nodes between J f and Jj _ That is, the sequence { Jj } is determined 
by availability (fig. 2). It IB = IE or MB = ME, then only one major axis (direction) exists.) 

the k lh designated node on the minor axis that intersects the major axis at Jj , where M k may not be 
monotonic or bounded for any i. The sequence {M k } is designated by the assignment of nodes in 
a cluster, if possible (fig. 2). 

the number of nodes ( 1 < Iq < N ; ) on the minor axis associated with Jj. The M j denotes the inter- 
section node of the minor axis with the major axis at Jj . 

the discrete function that first interrogates the availability of a path belonging to either RM h or RM V , 
which is initiated at the node (Jj, M^.), in the major direction. If applicable, the function then sub- 
sequently interrogates the availability of a path in the minor direction if, and only if, Jj + ( = IE or 
ME. The values of S ( J ; , M. ) are defined as follows: 

1 K i 

• NA = No availability (not any space on the major axis is free). 

• PA = Partial availability on the major axis. 

• TA = Total availability (both axes are free to complete the trace). 
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Figure 2. A representation of a typical trace. 



NI no further interrogation, implying total success or total failure. 


Definitions of Operations: 


C D the exchange of definitions of C and D. 

C D the exchange of definitions of C and D cannot take place. 

E[S(b, c)] evaluate S(b, c). 


INTRODUCTION 

An automated solution to the printed circuit board problem that is both efficient and cost effective continues to 
be in great demand. Most of the present solutions 1 ^ are either inexpensive but not efficient and fast, or efficient 
and fast but very costly. In either case, no solution exists that describes its algorithm symbolically with any degree 
of eloquence, nor are any algorithms sufficiently fast enough to facilitate parts placement for optimization. Many 
solutions are proprietary, so not much is written or known about the actual algorithms upon which these solutions 
are based. 

The purpose of this paper is to accomplish three objectives: 

• To facilitate optimization of arrangement of parts (for example, resistors and capacitors) on printed circuit 
boards for more sophisticated and expensive systems. 

• To make a very fast and efficient algorithm accessible at no cost for individual users or small companies. 

• To serve as a pedagogical tool that has heuristic value by presenting a formal symbolic structure that 
describes the entire algorithm. 

FORMAL SYMBOLIC STATEMENTS 

The purpose of the following algorithm, presented as a list of symbolic statements, is to construct a path from 

the source node, (IB, MB), to the target node, (IE, ME). A typical path (fig. 2) is described by the sequence of nodes 

(Jj, M 1( ), (J 2 , Mj 2 ), (J 2 , M 52 ), and (J 3 , M ]3 ). In this case, M 1; = M] 2 and M 52 = M] r In general, if (Jj, M k .) 

is an arbitrary element in a sequence leading to a solution, then the next element in the sequence will be either 

(J; , ,, M, ) or (J-, M/r, , ), where u > 0. Thus, the algorithm is as follows: 

1 + 1 *(i+ 1) 1 + u t| 

(a) Nj > 1, kj < Nj, and S(J jt M k ) = NA => E[S(J j( M (k + 1} )] 

(b) i > 1, Nj > 1, kj = Nj, andS(Jj,M k ) = NA => E[S(Jj_ p M (k+ ^ )] 

(c) i = 1, Nj > 1, kj = Nj, andS(Jj,M k ) = NA => (3 or (IB, MB) '5=^. (IE, ME) and E[S(J,, M , )] 

(d) i = 1, Nj > 1, kj = Nj,S(Jj,M k ) = NA,p^: p and (IB, MB) ^ (IE, ME) => NI 

(e) S(Jj,M k ) = PA=>E[S(Jj j,Mj )] 

(0 S(Jj,M k ) = TA=>NI 
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The efficiency of the algorithm is also predicated upon the process that determines the availability of rows and 
columns and determines the maintenance of the two respective directories. That is, a directory of rows, RM h , and 
a directory of columns, RM V , are maintained where each directory theoretically addresses different sides of the 
routing board. If required, any subsequent layers or surface mounts are handled similarly. 

In general, the determination of the availability at an arbitrary node, (p, 1), with respect to either RM h or RM V 
may be addressed accurately, keeping in mind that H p , and V p | are interchangeable as they relate to the following 
explanation (fig. 1). For the sake of specificity, H p f will be employed here. Thus, let (u Q , 1 0 ) be a target with re- 
spect to (p Q , 1 0 ) and where u Q > p 0 . Then the availability at (p 0 , 1 Q ) = min[H p j , u Q - p Q ] . However, if u Q < p 0 , 
then construct the equation T 0 = H , - (p Q - u 0 ) . Consider the two cases T Q > 0 and T 0 < 0. If T 0 > 0, then the 
availability at (p 0 , 1 Q ) with respect to (u Q , 1 0 ) is p 0 - u 0 . If T 0 < 0 , let T = H u { - (p 0 - u), where u Q < u < p Q . 
Then the availability at the node (p 0 , 1 0 ) that belongs to RM h is calculated using the bisection method whose di- 
rection of seek is governed by the sign of T. Moreover, because the maximum availability is sought, the number of 
iterations will be n = [LOG 2 (p 0 -u 0 )] + 1. 

The integrity of the algorithm, which is implemented by recursive descent, is not disturbed by either the defini- 
tion of the major and minor axes or the specificity of the terminal and initial nodes. Therefore, all four combinations 
of axes and nodes are transparent to the algorithm itself. Furthermore, an additional criterion to interchange axes or 
nodes could be predicated on the value of i or the magnitude of k ; for some i. 

The number of vias can be minimized from the contents of the directories (RM V and RM h ) after the solution 
is complete. This reduction results in a lower cost for the actual physical implementation of the solution. Finally, the 
updating of the directories dynamically is trivial, as the RM y and RM h contents are corrected based on the values 
of J- and the corresponding M k resulting from the monotonic sequence { J- 

The application of the theory presented confirms the objectives sought. Appendix A presents examples of typical 
solutions found. Moreover, a comparison study of the state-of-the-art algorithms is made. Appendix B addresses the 
efficiency, the cost, the ratio of solution to the Manhattan distance, and the central processing unit time. 

CONCLUSION 

An automated solution to the printed circuit board problem that circumvents the devastation of exponential growth 
as a function of complexity has been presented. Moreover, application of the theory displayed in Appendix A 
demonstrates its versatility and range. The statistics presented in Appendix B confirm the objectives sought. The com- 
puter program based upon the ideas presented here will be offered by COSMIC, located in Athens, Georgia. 

The author wishes to acknowledge Glenn Bever, Rodney Bogue, and Harry Chiles for their many helpful 
suggestions. 
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APPENDIX A 

EXAMPLES OF TYPICAL SOLUTIONS* 









































APPENDIX B 

TABLE OF COMPARISONS* 


Ratio, 


Algorithm 

actual distance/ 
Manhattan distance 

Completion 
rate, percent 

Cost, 

dollars 

CPU time, 
minutes 

Shape tech 

1.05 

98-100 

12,000 

10 

Pcad 

1.10 

92 

4,000 

30 

Hedgley 

0.97 

93 

0 

0.25 


Chart statistics are reflective of a one-layer board that is two-sided with equivalent grid sizes. 
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